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Qubit coherence decay down to threshold: influence of substrate dimensions 
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Keeping single-qubit quantum coherence above some threshold value not far below unity is a 
prerequisite for fault-tolerant quantum error correction (QEC). We study the initial dephasing 
of solid-state qubits in the independent-boson model, which describes well recent experiments on 
quantum dot (QD) excitons both in bulk and in substrates of reduced geometry such as nanotubes. 
Using explicit expressions for the exact coherence dynamics, a minimal QEC rate is identified in 
terms of the error threshold, temperature, and qubit-environment coupling strength. This allows us 
to systematically study the benefit of a current trend towards substrates with reduced dimensions. 

PACS numbers: 03.65.Yz, 78.67.Hc, 63.20.kd 
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Two-level systems with long-lived quantum coherence 
are candidate qubits, the units of quantum information 
[l| . Interactions between solid-state qubits and substrate 
phonons cause decoherence. It seems natural to fight de- 
coherence by reducing the number of substrate degrees 
of freedom, by going from bulk to planar or linear ge- 
ometries. In fact, QDs can nowadays be embedded in 
confined structures such as freestanding semiconductor 
membranes [2|,y| or nanotubes and nanowires [J,la|- Such 
substrates may allow tailoring of the phonon spectrum 
and, thus, controlling the qubit dephasing. However, it is 
not obvious whether fewer substrate dimensions do mean 
less decoherence. It is interesting to compare photolumi- 
nescence measurements of single QDs in bulk environ- 
ment [a, LZI, Is Isl and in nanotubes [10|, luj . In both cases, 
pure dephasing due to deformation-potential coupling to 
acoustic phonons is the dominant decoherence mecha- 
nism |12l . |l3[ since relaxation occurs on a much longer 
time scale. Using a Markovian master equation, i.e. ap- 
proximating the dephasing as exponential decay with a 
coherence time T2, one finds T2 = 00 for bulk and a fi- 
nite T2 time for ID substrates (details below); this would 
be an argument against reducing substrate dimensions. 
However, both in bulk and in reduced geometries, fast 
dephasing at short times has been observed as a broad 
background in spectra [6|, [tI, [sj, |9|, [ill, [ij]. According 
to 11| , it is this non-exponential decay that may hamper 
applications for quantum information processing (QIP) 
with ID substrates. But are these really less ideal? 

Future QIP devices will require built-in QEC, as 
some decoherence is inevitable. From the information- 
theoretical side come stringent requirements for fault- 
tolerant QEC: gate error levels e should be less than 10"'^ 
[15J . a value that may be relaxed in the future. Usually 
assumptions such as local and Markovian decoherence go 
into the derivation of e. One may criticize these [16|, llTI ] 
and go beyond them 18|. We start from the other end, 



with a given e and a realistic qubit-bath model. 

In this paper, we study the minimal rate Wqcc at which 
single-qubit errors should be corrected, requiring the co- 
herence to stay above the threshold value 1 — e. The 
beautiful experiments on QDs in bulk (3D) [g, [sj, ll9| 
and on nanotube excitons (ID) llll| are both well ex- 
plained by the so-called independent-boson model, and 
we employ its generic version. We present analytical ex- 
pressions for Wqoc and for the exact coherence dynamics, 
with the substrate dimension as a free parameter, that 
reproduce the measured rich dynamics and lead to new 
predictions: how can parameters best be changed so that 
error correction is needed less frequently. 

Model. — To account for the coupling of the QD to 
quantized lattice vibrations we follow Refs. [3, |20| and 
employ the independent-boson model with Hamiltonian 
H = Ho + Hqh where Hq = Aaz/2 + J^k ^k^k^k de- 
scribes the uncoupled system of QD and bosons, with 
Pauli matrix cr^ for the qubit and creation and annihila- 
tion operators 6j^ for phonons of mode k. We focus on 
the dominant deformation-potential coupling to acous- 
tic phonons with dispersion uj^ = Vs\k.\ where Vs denotes 
the sound velocity. The qubit-boson interaction iJqb = 
^'^z X]k5k(fek + feLk) is characterized by microscopic cou- 
plings (7k Of Dg y^\k\/y/2hpVvs in terms of a deformation 
potential Dg for the excitons, the mass density p and 
the volume V [7|, [iJ, |2l|. The more specialized model 
in Refs. [7|,lll| takes different deformation potentials and 
confinement lengths for electrons and holes into account. 
As shown below, our simpler model is accurate and has 
the major advantage that explicit expressions for the co- 
herence decay can be obtained. The phonon bath gives 
rise to the spectral density J(w) — J2k l5kP<^('^ — i's|k|) 
which for deformation-potential coupling and spherical 
QD symmetry with confinement length Ig has the form 
Js(aj) = asW^o;,^"^ exp(— w/wc) with cutoff frequency 
'j-'c = Vs/lg and dimensionless coupling strength as- 
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FIG. 1: (color online). Time evolution of coherence for 
bulk structures with s = 3 at different scaled temperatures 
= k-sT/hwc with az = 0.1. Inset: Non-monotonous 
temperature-dependence of the decay time f* for different 
coupling strengths as = O.Of (solid), 0.1 (dashed), and 0.8 
(dotted). i« is defined as the time at which the coherence has 
dropped 'half-way' to [1 -I- C3(oo)]/2, see Ref. [19| . 



Exact dynamics for arbitrary dimension. — We find the 
explicit expression for the qubit coherence Cs{t)^ i.e. the 
exciton polarization after a J-like excitation pulse. In 
terms of the scaled temperature 9 — k^T /HlJc, the dy- 
namics is given by Cs{t) = exp[— As(t)], with exponent 



A,(t) = 8a,(- 



F^'~^\e) - ReF^'-^\e[l -f iu^t]) 



/ cos[(s-l)arctan(^c^)] _ 
^'^ '^ (l + c.2t2)(-i)/2 



(1) 



where T{z) denotes Euler's Gamma function, F{z) = 
logr(0) and its n-th derivative F^"^ which for positive 
n equals the Polygamma function '^n-i{z). Hereby we 
generalize previous results for qubit dephasing [7|, |l7| to 
spectral densities with arbitrary s = 1, 2, 3, . . . 

We will now briefly discuss the qualitatively different 
dynamics for 3D, 2D, and ID substrates and present new 
analytical results, all based on Eq. ([1]), before comparing 
the error correction rates Wqec for the three geometries. 

3D substrates. — Dephasing due to acoustic phonons 
in bulk geometries is described by Eq. ([1]) for s = 3 (for 
deformation-potential coupling, the parameter s equals 
the dimension of the substrate). While a Markovian mas- 
ter equation for our model would predict no decoherence 
at all, the exact dynamics in Fig. [1] exhibits a fast initial 
decay of coherence, after which the coherence stabilizes to 
the final value 03(00) = exp[-8a3(6'^*i(6') - 1/2)], with 
zero-temperature limit exp(— 4q!3). Hence considerable 
initial decay may occur even at low temperatures, as the 
solid (blue) curve in Fig.[T]illustrates. Using the parame- 
ters of Refs. [TI, |9| for GaAs/InAs self-assembled QDs we 
obtain a typical value 03 = 0.8±0.3andwc = S-lO^^g-i^ 
in close agreement with the experimental curves [a, ll9| . 
For the exciton in the QD we find Z3 = 10 nm, a value in 
between the known confinement lengths of electron and 



hole wave functions. In Ref. 13| a decay time i, was 
defined as the time at which half of the coherence that 
finally will get lost is lost. Theory and experiment agreed 
well on the point that i* behaves non-monotonously as 
a function of temperature. In the inset of Fig. [1] we find 
for our model similar non-monotonous behavior, another 
illustration that dephasing of real QDs is also well de- 
scribed by our generic model. Notice the short time scale 
of the decay, i* < lj^^. There is monotonous depen- 
dence on other parameters: i* gets longer by making the 
QD larger, phonon velocity Vs slower, or the coupling a^ 
smaller. As Fig. [1] shows, the amount of coherence lost 
during short times is considerable for this realistic choice 
of parameters. If for QIP applications any phase errors 
beyond the percent level have to be corrected, then this 
should occur well within t, . 

Experimentally the coherence dynamics is usually in- 
ferred from absorption spectra. The Fourier transform 
of the highly non-exponential coherence dynamics in 
Fig. [1] predicts a highly non-Lorentzian spectrum. Since 
£3(00) > 0, the zero-phonon line at frequency A/S is even 
a delta function, with weight c^{oo) 12]. In practice this 
spectral line has a very narrow finite width, finite due to 
slow processes not described by the independent-boson 
model. However, with QIP in mind, we are more in- 
terested in the fast initial decay, which shows up as the 
broad background of the zero-phonon line. 

2D substrates. — Let us briefly also consider dephasing 
in planar geometries. To our knowledge, experiments 
analogous to ll|, ll9[ with planar geometries have not 
yet been carried out. The exact dynamics is obtained 
by setting s = 2 in Eq. H]). Again a fast initial decay 
is found, on a time scale u}~^ . The coherence C2(t) does 
not stabilize to a finite value, but vanishes algebraically 
ex l/(ciJci)^"^^j leading to a nontrivially broadened zero- 
phonon line in the absorption spectrum. 

ID substrates. — For s = 1, the spectral density 
Js{uj) becomes ohmic, i.e. linear in frequency up till the 
cutoff frequency Wc- This describes dephasing due to 
acoustic phonons in ID geometries such as nanotubes 
and nanowires, for which it is well known that the 
coherence dynamics strongly differs from the 3D case 
[111 ll4J . Detailed experiments have been performed only 
recently [11|. From Eq. ^ we find the exact analytical 
expression Ai(i) = 8ai{logr(6') - loglr(6'[l -I- iujct])\ - 
J log(l -I- ^2^2)}. For long times t > max{w-\ /i/fesT} 
this gives ci(t) = /cexp(— i/T2), i.e. exponential decay 
with coherence time T2 = h/AiraikBT. This also gives 
the zero-phonon line a finite width [ll|, [l4 ) . Master equa- 
tions would give the same T2 time, but not the prefac- 
tor K{ai,9) = [27r6l2«-7r2(6')]4"i by which the initial 
non-exponential decay remains noticeable also at long 
times. (See also the interesting discussion in Ref. [23], 
where another n is found.) In the limit <C 1 we find 
K = (27r6')^"i . Dephasing cannot be reduced indefinitely 
by lowering the temperature. Rather, the duration of 
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FIG. 2: (color online). Exact coherence dynamics for a ID en- 
vironment with ohmic spectral density with coupling strength 
ai = 0.01, for several scaled temperatures 9 = kBT/huic- The 
time-axis is scaled by 8. The dotted line shows the approxi- 
mated exponential decay as given by the T2 time (see text). 
Inset: same as main figure, but on a shorter time scale. The 
arrows indicate an error threshold e = 10"^. 



the non-exponential decay is increased. For the experi- 
ments in Ref. , 1 Ij with isolated single- wall nanotubes we 
estimate ai = 0.1 ± 0.05 and cj^ = 20 • lO^^g-i -y^jth 
temperatures ranging from 5K to 32K, i.e. 9 between 
0.033 and 0.21, k assumes values between 0.49 and 0.92. 
Thus non-exponential dephasing is important in state-of- 
the-art ID systems, the more so for lower temperatures. 
Again, our main interest is the initial decay itself, 
rather than its effect at long times. For < 1, the 
time scale of the non-exponential decay is h/k^T, in clear 
contrast to the temperature-independent t* for bulk sys- 
tems. Figure [2] shows the typical coherence dynamics of 
QDs on ID substrates. We chose as = 0.01, an order 



of magnitude smaller than in the experiment ll|, but 
non-exponential decay would be important even then. 
The figure shows curves for three different tempera- 
tures, as well as their master-equation approximations 
exp(— t/r2). The latter coincide due to scaling of the 
time axis. By contrast, the three exact curves do not co- 
incide at all: the low-temperature curve (9 = 0.01) is sys- 
tematically lower and the high-temperature curve higher 
than their exponential-decay approximations. The inset 
of Fig. [2] shows that when asked how long ci{t) manages 
to stay above 0.999, the exponential-decay curve may be 
too optimistic, too pessimistic, or accurate by chance. 

Rate of QEC. — There is the danger of comparing ap- 
ples and oranges when studying the effect of substrate 
dimensions. Fortunately, in state-of-the-art experiments 
both as and ai turn out to be of order 10"^ (see above). 
For current and future experiments, it is therefore useful 
to depict the as-dependence of dephasing for all three 
geometries in the same range as — 0.1 and smaller. 

Let us now assume that phase errors have to be cor- 
rected if coherence drops from Cs(0) = 1 to Cs (27r/wqoc) = 
(1 — e) for some error threshold e, i.e. from the exact de- 
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FIG. 3: (color online). Temperature dependence of the quan- 
tum error correction rate u^cc for substrate geometries with 
spectral densities Js{(j->) with different s but equal couplings 
Qs, for several threshold values e. In the main figure as = 0.1, 
and as = 0.001 in the inset, where case s = 2 is not shown. 
The y-axes are scaled by yfl. Curves are based on Eq. ([l]). 



phasing dynamics ([T]) we identify the minimal rate Wqoc 
at which phase errors have to be corrected in order to 
preserve coherence in an idle qubit. To be optimistic, 
we assume that in each step the error can be corrected 
perfectly and instantaneously. Central idea is that struc- 



tures with lower rates 



^qcc 



are better suited for the im- 



plementation of QEC. In Figure [3] we compare (jjqec as a 
function of temperature for bulk, planar, and linear sub- 
strate geometries. All couplings are as = 0.1. Notice 
that Wqcc on the vertical axis is scaled by ^/e. We find 
that the curves for e — 10^^, 10^'^, and 10"^ overlap for 
all three geometries. Thus Wqoc scales as y^, at least for 
current experimental couplings as = 0.1 and e < 10^^. 
However, this simple scaling breaks down for very weak 
coupling as = 0.001, as shown in the inset of Fig. [3l 
Scaling with y^ only holds if the parabolic short-time 
approximation Cs{t) « 1 — r]st^ is still valid at the time 
that Cs{t) assumes the threshold value (1 — e). We find 
ris{as,e) = 2as[2{-ey+'^'^s{0) - si], and for smaU e the 
error correction rate Wqoc = 2Tr^r]s/e. For as = 0.1, 
the temperature dependence of rjs and hence of Wqec is 
completely negligible, even up to temperatures as high 
as hujc/kB. This corroborates that the first stage of pure 



dephasing is due to vacuum noise [Hl, |23|. Fig. [3] also 
shows that the rates for ID are smaller than for 3D or 2D 
geometries, although by factors less than 10. Generally 
we find the central result that linear substrate geometries 
like nanotubes and nanowires will perform best. All rates 
<^qcc in Fig- 131 a-re high, in between cutoff lUc and qubit 
frequency A/h ~ 10^^ s"^, since already the fast initial 
decay reaches the error threshold. 

It would be very challenging to implement a quan- 
tum error correction protocol for phase errors at such 
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FIG. 4: (color online). Quantum error correction rate ujqcc 
as a function of coupling Us, for several error thresholds e. 
Curves were made using Eq. ([T}. For all three pairs of curves, 
s = 1 corresponds to the lower, and s = 3 to the upper 
curve. The curves for s = 3 stop at critical coupling strengths 
indicated by dots: For smaller 03, the coherence C3(t) never 
decays below threshold. The temperature is T = O.Olhcjc/kB- 



high rates. In our model, lower rates Wqcc could be re- 
ahzed by reducing the cutoff frequency tOc or the cou- 
phngs as ■ Figure [3] shows the effect of the latter strat- 
egy. For the largest couplings as — 0.1, the message is 
as in Fig. [3l given an error level e, linear substrates have 
lower rates LUqec than planar or bulk substrates. As the 
as are decreased, this message is essentially unchanged, 
until suddenly for bulk substrates a^ gets so small that 
03(00) > 1 — e, i.e. the final coherence stabilizes above 
the error threshold. In that situation - which can not 
occur for ID substrates - Wqoc vanishes: no error cor- 
rection is needed. State-of-the-art exciton qubits in QDs 
have couplings that are more than one order of magni- 
tude larger than the largest critical coupling shown {i.e. 
for e = 0.01). Obviously Wqcc would also become smaller 
if larger errors e were allowed. The challenge is here to 
come up with QEC protocols that tolerate larger faults. 
All in all. Fig. |4] shows that for fixed as = a and Wc, 
linear substrate geometries are to be preferred for their 
lower ( 



^qec? 



unless couplings can be substantially reduced. 

Discussion and conclusions. — Inspired by recent mea- 

we have analyzed the first few per- 



UM 



surements 

cents of loss of quantum coherence of a solid-state exci- 
ton qubit on 3D, 2D, and ID substrates. It is mainly 
this initial decoherence that is important for QIP appli- 
cations when supplemented with fault-tolerant QEC. We 
proposed and focused on the important quantity Wqcc, 
the minimal rate at which quantum errors should be cor- 
rected. Its temperature dependence turned out to be 
negligible. For QD exciton qubits in a bulk substrate, 
the coherence may stabilize above the threshold, as for 
spin qubits [2J] , but corresponding couplings as are cur- 
rently not weak enough. 

Let us return to our initial question: Is it beneficial 
for QIP applications to reduce the dimensions of the 



substrate, as far as dephasing is concerned? From a 
master-equation perspective it is not, and worries about 
non-exponential decay were expressed especially for ID 



structures [11|. We have presented analytical solutions 
of a generic but accurate independent-boson model to 
show how fast initial dephasing occurs for ID, 2D, and 
3D geometries alike, and identified the minimal rate at 
which single-qubit error have to be corrected. In all cases 
the rates are high, in between the cutoff and the qubit 
frequency, which poses a challenge for QIP applications. 
However, we found that qubits on ID substrates require 
the lowest error correction rates. 

This work has been supported by the the DFG through 
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